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1. Executive Summary 

We present a set of compelling science cases for the 31-45 GHz Band 1 ALMA Receiver. 
The science ranges from nearby stars to the re-ionization edge of the Universe. Two of the 
science cases provide additional leverage on the present ALMA Level One Science Goals and 
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are seen as particularly powerful motivations for building the Band 1 Receiver: (1) detailing 

the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, 
requires continuum observations out to ~ 30 GHz (~ 1 cm); (2) detecting CO 3-2 spectral 
line emission from Galaxies like the Milky Way during the era of re-ionization, 6.5 < ^ < 10 
also requires Band 1 Receiver coverage. 
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2. Introduction 

The Atacama Large Millimeter/submillimeter Array (ALMA) will be a single research 
instrument composed of at least fifty-four 12-m and twelve 7-m high-precision antennas, 
located at a very high altitude of 5000 m on the Chajnantor plain of the Chilean Andes. 
The weather conditions at the ALMA site will allow transformational research into the 
physics of the cold Universe across a wide range of wavelengths, from radio to submillimeter. 
Thus, ALMA will be capable of probing the first stars and galaxies and directly imaging 
the disks in which planets are formed. ALMA will be a complete astronomical imaging and 
spectroscopic instrument for the millimeter/submillimeter regime, providing scientists with 
capabilities and wavelength coverage that complement those of other research facilities of its 
era, such as the Expanded Very Large Array (EVLA), James Webb Space Telescope (JWST), 
Thirty Meter Telescope (TMT), European Extremely Large Telescope (E-ELT), and Square 
Kilometer Array (SKA). ALMA will be the preeminent platform for astronomical research 
at millimeter (mm) and submillimeter (submm) wavelengths for decades to come. 

ALMA will revolutionize many areas of astronomy and an amazing breadth of science 
has already been proposed (see for example the ALMA Design Reference Science Plan). The 
technical requirements of the ALMA Project are, however, driven by three specific Level 
One Science Goals: 

(1) The ability to detect spectral line emission from CO or C-|- in a normal galaxy like the 
Milky Way at a redshift oi z — 3, in less than 24 hours of observation. 

(2) The ability to image the gas kinematics in a solar-mass protostellar/ protoplanetary disk 
at a distance of 150 pc (roughly, the distance of the star-forming clouds in Ophiuchus or 
Corona Australis), enabling one to study the physical, chemical, and magnetic field structure 
of the disk and to detect the tidal gaps created by planets undergoing formation. 

(3) The abihty to provide precise images at an angular resolution of 0.1 arcseconds. Here 
the term "precise image" means an accurate representation of the sky brightness at all points 
where the brightness is greater than 0.1% of the peak image brightness. This requirement 
applies to all sources visible to ALMA that transit at an elevation greater than 20°. 

ALMA was originally envisioned to provide access to all frequencies between 31 and 
950 GHz accessible from the ground. During the re-basehning exercise undertaken in 2001, 
the entire project was scrutinized to find necessary cost savings and the two lowest receiver 
frequencies, covering 31-84 GHz, were among those items delayed beyond the start of science 
operations. Nevertheless, Band 1, covering 31-45 GHz, was re-affirmed as being a high 
priority future item for ALMA. 

In May 2001, John Richer and Geoff Blake prepared a document Science with Band 1 
(3I-45 GHz) on ALMA as part of the re-baselining exercise. Key arguments for maintaining 
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the receiver were put forward and included: (1) exciting science opportunities, bringing in a 
wider community of users; (2) a significantly faster imaging and survey instrument than the 
upgraded EVLA, especially due to the larger primary beam; (3) access to the southern sky at 
these wavelengths; (4) excellent science possible even in "poor" weather; (5) a relatively cheap 
and reliable receiver to build and maintain; and (6) a very useful engineer ing/debugging tool 
for the entire array. 

The Richer/Blake document was followed by an ASAC Committee Report in October 
2001 in which the addition of Japan into the ALMA project re-opened the question of 
observing frequency priorities for those receiver bands which had been put on hold during 
re-baselining. The unanimous recommendation of the ASAC was to put Band 10 as top 
priority, followed by a very high priority Band 1. The key science cases for Band 1, at that 
time, were seen to be (1) high resolution Sunyaev-Zel'dovich (SZ) imaging of cluster gas at 
all redshifts; and (2) mapping the cold ISM in Galaxies at intermediate and high redshift. 

The scientific landscape has changed significantly since 2001 and thus it is time to update 
the main science drivers for Band 1. The ALMA Development Plan also now lies on the 
horizon, and those parties interested in longer wavelength observing with ALMA recognize 
that now is the time to put forward their best cases. Thus, in the Fall of 2008, two dozen 
astronomers from around the globe met in Victoria, B.C. to discuss this issue. This paper 
represents both a record of that meeting and an attempt to move toward concrete proposals 
for outstanding science with Band 1. In Section [3] we present two science cases that reaffirm 
and enhance the already established ALMA Project Level One Science Goals. In Section 
|4] we provide a selection of science cases that reinforce the breadth and versatility of the 
Band 1 Receiver. Section [5] discusses both weather considerations at the ALMA site and 
the complementarity of ALMA versus the EVLA at Band 1 frequencies. Finally, Section |6] 
summarizes the paper. 



3. Level One Science Cases for Band 1 



In this section we present two science cases that reaffirm and enhance the already 
established ALMA Project Level One Science Goals: Evolution of Grains in Disks Around 
Stars (3.1) and The First Generation of Galaxies (3.2). 



-5 - 



3.1. Evolution of Grains in Disks Around Stars 

3.1.1. Protoplanetary Disks 

Planet formation takes place in disks of dust and gas surrounding young stars. It 
is within these gas-rich protoplanetary disks that dust grains must agglomerate from the 
sub-micron sizes associated with the interstellar medium to larger pebbles, rocks and plan- 

etesimals, if planets arc ultimately to be formed. The timescale of this agglomeration process 
is thought to be a few tens of Myr for terrestrial planets, while the process leading to the 
formation of giant planet cores remains uncertain. Core accretion models require at least 
a few Myr to form Jovian planets (Pollack et al. 1996), while dynamical instability models 
could form giant planets on orbital timescales (t ^ 1 Myr; Boss 2005). 

Gravitational instability models require high disk masses in order to form planets. So 
far, most accurate disk mass estimates come from submm and mm observations, where the 
dust is optically thin. Andrews & Williams (2007a, 2007b) show that submm observations of 
dozens of protoplanetary disks reveal that only one system could be gravitationally unstable, 
conflicting with the high frequency of Jovian planets seen around low mass stars. Have these 
relatively young (1-6 Myr) systems already formed planets, or is most of the dust mass 
locked into larger grains and therefore not accounted for in submm and mm observations? 
If grain growth to centimeter sizes has occurred, most of a disk's dust mass would reside in 
the large particle population, which would emit at longer mm and cm wavelengths. Greaves 
et al. (in prep.) compare disk masses for objects in Taurus and Ophiuchus derived from 3 
and/or 7 mm with 1.3 cm dust fluxes. The cm-based masses are found to be higher, by 
factors up to 5, than the shorter mm-based values, indicating that a significant fraction of 
the disks' total dust masses are indeed locked into cm-sized grains. 

Identifying where and when dust coagulation occurs are critical to constrain current 
models of planetary formation. Growth from sub-micron to micron-sized particles can be 
traced with infrared spectroscopy and imaging polarimetry. The next step, growth beyond 
micron sizes, is readily studied by determining the slope of the spectral energy distribution 
(SED) of the dust thermal emission at submm and mm wavelengths. The dust mass opacity 
index at wavelengths longer than 0.1 mm is approximately a power-law whose normalization 
depends on the dust properties, such as composition, size distribution, and geometry (Draine 
2006). The index of the power law is (3. The presence of large grains is detectable through a 
decrease in the index /3, which can be derived directly from the slope of the Rayleigh- Jeans 
tail of the SED, a, where /3 = a — 2, when the emission is optically thin. Early studies of 
disks revealed that the (3 values of disks are substantially lower than the typical ISM value 
of ~ 2 (e.g., Weintraub et al. 1989; Adams et al. 1990; Beckwith et al. 1990; Beckwith & 
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Sargent 1991, Mannings & Emerson 1994). 

The key stumbling block to the interpretation of /3 occurs when the disk is not resolved 
spatially. The amount of flux detected at a given wavelength is a function of both /3 and the 
size of the disk (Testi et al. 2001). Resolving the ambiguity therefore is truly a matter of 
resolution, and sufficient resolution is only offered at these wavelengths by interferometers. 

Among the three high level science goals of ALMA is the ability to detect and image gas 
kinematics in protoplanetary disks undergoing planetary formation at 150 pc. At ALMA's 
observing wavelengths, its capability for imaging the continuum dust emission in these disks 
is also second-to- none. The re-baselined ALMA, however, will only reach a wavelength of 3.6 
mm. Given that dust particles emit very inefficiently at wavelengths longer than their sizes, 
the re-baselined ALMA design will not be sensitive to particles larger than ~ 3 mm. This 
negates ALMA 's potential ability to follow the dust grain growth from mm-sized to cm-sized 
pebbles in protoplanetary disks. 

Figure [T] shows the SEDs for three different circumstellar disk models, computed using 
the full dust radiative transfer MCFOST code (Pinte et al. 2006; Pinte et al. 2009 ). The 
model parameters are representative of protoplanetary disks (although there is substantial 
object-to-object variation). The circumstellar disk is passively heated by a 4000 K, 2 Lq 
central star and the system is located 160 pc away. The dust component of the disk is 
assumed to be fully mixed with the gas and the latter is assumed to be in vertical hydrostatic 
equilibrium. The disk extends radially from 1 AU to 100 AU. The total dust mass in the 
model is 10"'^ M© (the gas component is irrelevant for continuum emission calculations, so 
its mass is not set in the model, though a typical 100:1 gas:dust ratio is generally assumed). 
The dust population is described by a single power law size distribution N{a) oc a~^-^ with 
a minimum grain size of 0.03 /zm and extending to 10 /im, 1 mm or 1 cm depending on the 
model. The dust composition is taken to be the "astronomical silicates" from Draine (2003). 

Figure [T] reveals that observations in the ALMA Band 1 spectral region are crucial for 
determining whether grain-growth to cm-sizes is indeed occurring. The 1 cm flux density of 
the max^j^e =1 cm disk model is ~ 50 /xJy, comparable to the 1-sigma sensitivities provided by 
ALMA's Band 1 with 1 minute integration. Besides ALMA, there are no existing or planned 
southern astronomical facilities capable of observing to such depths at these frequencies. 
Therefore, if Band 1 is not built there will be no way of putting ALMA observations of 
protoplanetary disks in the context of coagulation of dust grains to centimeter sizes. 

By complementing observations in other ALMA Bands, Band 1 will provide a crucial 
longer wavelength lever to minimize the uncertainty in a. Evidence for small pebbles has 
been detected in several disks (Rodmann et al. 2006). The prime example is TW Hya, a 
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protoplanetary disk 50 pc from the Sun (Wilner et al. 2000). Its SED is well matched by 
an irradiated accretion disk model fit from 10s of AU to an outer radius of 200 AU and 
requires the presence of particle sizes up to 1 cm in the disk (see Figure |2|. The measured 
P is 0.7 ± 0.1 (Calvet et al. 2002). To date, no trend in /3 has been detected with stellar 
luminosity, mass or age. Lower a values are associated with less 60 /xm excess, however, 
suggesting that settling or agglomeration processes could be removing the smallest grains, 
decreasing the shorter wavelength emission (Acke et al. 2004). 

At the resolution provided by its longest baselines at 31 GHz (~ 0.1 arcseconds), ALMA 
will easily resolve protoplanetary disks at the distance of the closest star-forming regions 
(50-150 pc). These resolved images will provide the most accurate determination of the 
disk's dust mass. The dust distribution at cm wavelengths can then be compared to mm 
and submm images, revealing where in the disk dust coagulation is occurring. Through 
detection of concentrations of large grains, protoplanets in formation can be identified. Such 
condensations are predicted by simulations of gravitational instability models (see Figure |3^; 
Greaves et al. 2008) and have been detected in the nearby star HL Tau (Figure |3]d; Greaves 
et al. 2008). ALMA's high resolution is clearly critical to the detection of such substructures 
within protoplanetary disks. 

Detecting dust emission at cm wavelengths also requires high sensitivity, because its 
brightness is several orders of magnitude lower than in the submm. In addition, at wave- 
lengths longer than 7 mm, the contribution from other radiative processes, such as ionized 
winds, can contribute significantly to the total flux and complicate the interpretation of 
detected emission. Rodmann et al. (2006) found that the contribution of free-free emission 
to the total flux is typically 25% at a wavelength of 7 mm. The detailed sampling of the 
31-45 GHz (1 cm-6 mm) spectral range enabled by Band 1 would increase substantially the 
sampling rate in the region where emission is detected from both the free-free and thermal 
dust emission components. The synergy with the EVLA will provide a longer wavelength 
lever for sources observed in common, providing an estimate of the free-free contribution 
to the Band 1 flux. This would not be essential, however, since detailed sampling of these 
data are crucial to quantify accurately the relative amounts of free-free and continuum emis- 
sion, and thereby determine precisely the thermal emission due to large dust grains (i.e., 
protoplanetary material). 



3.1.2. Debris Disks 

A different process is at work in the production of pebble-sized bodies in evolved, gas- 
poor disks around main-sequence stars. Destructive, coUisional cascades from even larger 
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planetesimals through to cm, mm and then micron-sized particles provides ongoing replen- 
ishment of the debris population (Wyatt 2009; DuUemond & Dominik 2005). The methods 
of detection of large (i.e., cm-sized) grains is the same as in protoplanetary disks, despite 
their origin in destructive, rather than agglomerative, processes. In each case, the longer 
the wavelength at which continuum emission is detected, the larger the grains that must be 
present in the system. 

Band 1 images could show higher contrast features in debris disks compared to other 
ALMA bands, due to the longer resonant lifetimes of the larger particles that dominate 
the emission. This in turn will help detect any edges and gaps in the disks. Dramatic 
changes in the morphology of debris disks as a function of wavelength have already been 
observed (e.g., Maness et al. 2008), but not yet at such long wavelengths as Band 1 will 
probe. When observed, such structures are often considered signposts to the existence of 
planets. These detections will be challenging compared with detecting forming condensations 
in protoplanetary disks. This is because of the relatively low surface brightness and large 
spatial distribution of the debris disks, which can be 100s of AU in radii and can be found 
much closer to the Sun than protoplanetary disks. The closest disks could subtend as much 
as ~ 150 arcseconds on the sky (assuming a 300 AU diameter disk at 2 pc). Therefore, 
ALMA's larger field of view than other long wavelength instruments, such as the EVLA, will 
be very advantageous for imaging these objects, although mosaicking will still be required to 
image the largest ones on the sky. ACA observations will be instrumental in this research. 



3.2. The First Generation of Galcixies: 
Moleculcir gas in galetxies during the epoch of reionization 

The first generation of luminous objects in the Universe began the process of re-ionizing 
the intergalactic medium (IGM). The detection of large-scale polarization in the cosmic 
microwave background (CMB), caused by Thomson scattering of the CMB by the IGM 
during re-ionization, suggests that the Universe was significantly ionized as far back as 
z ~ 11. 0± 1.4 (Dunkley et al. 2009). The "near" edge of the era of re-ionization has been 
inferred from the detection of the Gunn-Peterson effect (Gunn & Peterson 1965) toward 
galaxies with -2 > 6 (Fan et al. 2006a,b). The nearly complete absorption of all continuum 
shortward of the Lya break is due to moderate amounts of neutral hydrogen in the IGM, 
suggesting re-ionization was complete by z ~ 6. The Gunn-Peterson effect also insures that, 
at these redshifts, the Universe is opaque at wavelengths shortward of ~ 1 /xm. 

To study the first generations of galaxies, and to understand the origins of the black 
hole-bulge mass relation, it will be necessary to study galaxies in the 6 < 2; < 11 range. 
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Quasar hosts and other sources are rapidly being discovered at the near end of this range 
(e.g. Cool et al. 2006; Mortlock et al. 2008; Glikman et al. 2008; Willott et al. 2009), 
and searches are underway for even more distant objects (e.g. Ota et al. 2008; Bouwens et 
al. 2009). 

Recently, CO has been detected in galaxies at redshifts > 6. These and other obser- 
vations in the cm/mm of 2; > 6 galaxies are summarized by Carilli et al. (2008). Current 
instrumentation sensitivities are such that detections are limited to hyperluminous infrared 
galaxies, i.e. Lpm > 10^^ Lq. Only a small fraction of galaxies are this luminous. The 
best-studied such object is J1148+5251 with a redshift of 2; = 6.419 (see Carilli et al. 2008). 
For example, Walter et al. (2004) imaged the CO J=3-2 emission (Figure |4]) using the VLA, 
from which they were able to infer the dynamical mass. Walter et al. (2009) were not able 
to detect the [Nil] 205^111 line, but did detect the CO J=6-5 transition. These and other (dust 
continuum) observations show that there was already a significant abundance of metals and 
dust by this epoch. 

Figure [5] shows the frequency of each transition of ^^CO, from J=l-0 through J=10-9, as 
a function of redshift. Also shown are the frequency ranges of the ALMA Bands (excluding 
Band 2). In the redshift range of the era of re- ionization (2; > 6), Band 3 receivers would be 
able to detect J=6-5 emission in the range 4.8 < 2; < 7.2. Band 1, however, is almost ideal, 
since its redshift coverage for J=3-2 emission is approximately 6.5 ^ 2; < 10. If the Band 1 
coverage were to be extended to 50 GHz, the redshift coverage would expand to 6 < 2; < 10. 

3.2.1. Quasar Host Galaxies 

The discovery of molecular gas in quasar host galaxies at z ~ 6, when the Universe 
was less than 1 Gyr old (Waher et al. 2003; Bertoldi et al. 2003; Carilli et al. 2007), has 
opened a new window on the study of gas in systems that contributed to the reionization 
of the Universe. Studies of how the molecular gas properties should evolve, and how they 
can be used to reveal the dynamics of these massive systems, have recently prompted a new 
generation of semi-analytic models with the further aim of understanding how high-redshift 
quasars fit within the context of large-scale structure formation. Li et al. (2007, 2008) have 
used state of the art N-body simulations to show that the observed optical properties of high- 
redshift quasars can be explained if these objects formed in the most massive dark matter 
haloes (~ 8 x 10^^ M©) early on. These models predict that the most luminous quasars 
should evolve due to an increase of major mergers, which one would expect to find evidence 
for in the CO line profiles and the spatial distribution of the molecular gas (Narayanan et 
al. 2008). Detailed radiative transfer models of the FIR spectral energy distribution of these 
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systems have been driven by the observations of one z — 6.42 quasar (namely J1148+5251; 
Walter et al. 2003, 2004), and so larger samples of CO-detected quasars are needed to provide 
better constraints on the models. Current 3 mm surveys of high-J CO line emission in 2; ~ 6 
FIR-luminous quasars are being conducted with the PdBI, having successfully detected CO 
line emission in eight objects (Wang et al. 2009). 

The availability of Band 1 receivers on ALMA would permit complementary surveys 
and imaging of CO J=2-l line emission in southern quasars out to z < 6.4, and CO J=3-2 
line emission for those at 2; > 6.7. Assuming a 150 /uJy CO J=2-l line of width ~600 km s~^ 
at 2; = 5.7, a 5(T detection would take less than 4 hours with 50 ALMA antennas. 

3.2.2. Lyrnan-a Emitters 

The rarity of the luminous quasars at early times suggests that their UV emission is 
unlikely to contribute significantly to the reionization of the Universe (e.g. Fan et al. 2001). 
A more important type of galaxy in the context of cosmic reionization, and a significant 
fraction of the star-forming galaxy population at 2; ~ 6 are the Lyman-a emitters (hereafter 
LAEs), which have been discovered through their excess emission in narrow-band images 
centered on the redshifted Lyman-a line (e.g. Hu et al. 1998; Rhoads et al. 2000; Taniguchi 
et al. 2005). While the star-formation rates in LAEs, inferred from their UV continuum 
emission, are a few tens of solar masses per year (e.g. Taniguchi et al. 2005), their number 
density and the shape of the Lyman-a emission line provide important probes of physical 
conditions in the Universe around the epoch of reionization. As such, it is of fundamental 
importance that we understand the properties related to their star-formation activity. In 
particular, we need to quantify the amount of molecular gas available for fuel. Wagg, Kanekar 
& Carilli (2009) have used the Green Bank Telescope to search for CO J=l-0 line emission 
in 2; > 6.5 LAEs, targeting two sources including the highest spectroscopically confirmed 
redshift LAE z — 6.96 (lye et al. 2006). The limits to the CO fine luminosity impfied 
by the non-detections of CO J=l-0 in these two objects suggest modest molecular gas 
masses (;SlO^° M©). This conclusion, however, is based on observations of only two objects, 
and future studies would benefit from the sensitivity gained by observing higher-order CO 
transitions, whose flux density may scale as z/^ due to a contribution to the molecular gas 
excitation by the cosmic microwave background radiation (19 K at 2; = 6). At these redshifts, 
such studies would require the sensitivity of the Band 1 receivers on ALMA. 
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4. A Broad Range of Science Cases 

Along with tlie two science cases presented above, tliere is a wealth of scientific oppor- 
tunity available to the ALMA community if the Band 1 Receiver is built. Here we highlight 
a selection of science cases which would significantly benefit from Band 1 on ALMA. 

4.1. Continuum Observations with ALMA Band 1 

The astrophysical continuum radiation at wavelengths of ~ 1 cm is relatively unexplored. 
Yet this radiation is key to understanding radio emission mechanisms and probing regions 
that are optically thick at shorter wavelengths. The resolution and sensitivity of ALMA Band 
1 will allow: (1) a diagnostic of the smallest interstellar dust grains, (2) an understanding of 
the nature of pulsar wind nebulae, (3) the detection of radio SNe, with constraints on stellar 
precursors and remnants, and (4) a diagnostic of X-ray binaries. 

4-1.1. Very Small Grains and Spinning Dust 

The last decade has seen the discovery of surprisingly bright cm-wavelength radio emis- 
sion in a number of distinct galactic objects but most notably in dark clouds (e.g., Finkbeiner 
et al. 2002; Casassus et al. 2008 (see Figure [6]); Scaife et al. 2009a). The spectrum of this new 
component of continuum radiation can be explained by electric dipole radiation from rapidly 
rotating ("spinning") very small dust grains (VSGs), as calculated by Draine & Lazarian 
(1998; DL98). It has also been seen as a large scale foreground in CMB maps, its spectrum 
peaking at ~40 GHz, spatially correlated with thermal dust emission. 

All of the existing work aimed at diagnosing this continuum emission is derived from 
CMB experiments on large angular scales, where the bulk of the radio signal occurs, and 
which will soon be surveyed by the Planck satellite. Details on small angular scales are 
crucial for probing star formation and circumstellar environments. Simply, progress in the 
understanding of the solid and gas states of the ISM requires sufficient resolution to separate 
the distinct environments. Directly measuring the VSG abundance and solid state physics 
is very exciting becasue VSGs play a central role in the chemical and thermal balance of 
the ISM. For example, the smallest grains account for most of the surface area available for 
catalysis of molecular formation. 

DL98 proposed that the grain size distribution in their spinning dust model would be 
dominated by VSGs, thought to be mostly PAH nanoparticles. The size distribution of VSGs 
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is poorly known as studies of interstellar extinction are relatively insensitive to its details. 
The existence of VSGs is supported by assertions that the existence of a significant amount 
of carbonaceous nanoparticles in the ISM could explain observations of the unidentified 
IR emission features and the strong mid-infrared emission component seen by IRAS which 
must result from starlight reprocessing of ultrasmall grains. The fraction of the ISM carbon 
content proposed to exist in VSGs considerably exceeds that fraction implied by the MRN 
dust size distribution, which is known to underestimate it. 

Obscrvationally determining PAH content in dust clouds is not straightforward. Where 
there is a strong source of UV flux present it is possible to identify PAHs by their spectral 
emission features. In the case of pre-stellar and Class clouds, however, these features are 
absent. With microwave observations from ALMA Band 1 constraining the spinning dust 
SED at similar resolution to, e.g., Spitzer or the forthcoming MIRI instrument on the JWST, 
it will be possible to measure the size distribution of VSGs directly from the data. This work 
will also be important in the context of circumstellar and protoplanetary disks, where the 
proposed population of VSGs may have important implications for disk evolution. Certainly, 
spinning dust emission will provide a better measure of the small grain population within 
circumstellar disks than PAH emission since favorable excitation conditions for PAHs exist 
only in the outermost layers of the disk. Since all the VSGs in the disk should contribute 
spinning dust emission, such emission will provide a much better probe of the mass in VSGs. 
Combining this information with the PAH emission features would then also give us a useful 
measure of sedimentation in disks. 

Theoretically, the spinning dust emission from a VSG population has been shown to 
dominate the thermal emission from stellar disks (HAeBe stars) at significant factors in 
the region < 50 GHz (Rafikov 2006). The existence of these VSGs has been confirmed 
obscrvationally using PAH spectral features as seen in the disks of Herbig Ae/Be stars (Acke 
& van den Ancker 2004) but it has not been detected in protoplanetary disks due to a lack 
of strong UV flux. Since spinning dust emission has been observed to be spatially correlated 
with PAH emission (Scaife et al. 2009b) spinning dust may provide a unique window on the 
small grain population of these disks. This is important in the context of disk evolution as 
these recent measurements conflict with the established view that dust grains are expected 
to grow as disks age. It may be the case that dust fragmentation is important in disks 
(DuUemond & Dominik 2005), or there is a separate population of very small carbonaceous 
grains, distinct from the MRN distribution (Leger & Puget 1984; Draine & Anderson 1985). 
This second proposition has not only important implications for the study of circumstellar 
disks but also for the complete characterization of dust and the ISM more generally. 

The arcsecond resolution necessary for these measurements will be achievable with a 
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number of the ALMA Band 1 configurations. From the models of Rafikov (2006), the 
difference between a thermal dust spectrum with /3 ^ 1 and the predicted spinning dust 
contribution for a brown dwarf disk would be observable at 5 a in a matter of minutes 
with ALMA Band 1. With longer observation times and consequently higher sensitivity, it 
will also be possible to distinguish between different grain size distributions and physical 
conditions within the disk (such as grain electric dipole moments, rotational kinematics, 
optical properties and catalysis of molecule formation) as well as gaseous physical conditions 
(gas temperature, density and ionization fraction). 

Spinning dust emission provides a unique insight into the VSG population under con- 
ditions where it is not possible to observe using mid-IR emission. The high resolution and 
excellent sensitivity of ALMA are ideal for differentiating the distinct environments where 
the VSG population resides and will be crucial for probing star formation and circumstellar 
regions. Band 1 will allow routine surveys of the new continuum component at its spectral 
maximum. 



4.1.2. Pulsar Wind Nebulae 

Pulsars generate magnetized particle winds that inflate an expanding bubble called a 
pulsar wind nebula (PWN) whose outer edge is confined by the slowly expanding supernova 
ejecta. Electrons and positrons are accelerated at the termination shock some 0.1 pc distant 
from the pulsar. Those relativistic particles interact with the magnetic field inside the wind 
blown bubble to produce synchrotron emission across the entire electromagnetic spectrum. 
Particles accelerated at the shock form toroidal structures, known as wisps, and some of 
them are collimated along the rotation axis of the pulsar, contributing to the formation 
of jet-like features. The synchrotron emission structure in the post-shock and jet regions 
provide direct insight on the particle acceleration process, magnetic coUimation, and the 
magnetization properties of the winds in PWNe. These observations have so far (except for 
the Crab Nebula) been limited to the X-ray band with the Chandra satellite (e.g. Helfand 
et al. 2001). 

ALMA has the sensitivity and resolution neccessary to detect PWNe features at high 
radio frequencies, where we can detect the emission from relativistic particles which have 
much longer lifetimes than in X-rays. At cm/mm- wavelengths, flat-spectrum synchrotron 
PWNe stand out over steep-spectrum SNRs [e.g. as seen in the Vela PWN, Hales et al. (2004) 
and also as discussed in Sec. 4.1.3 below, and illustrated in Fig. [7], taken from Bietenholz, 
et al. (2004)] with minimal confusion from the Rayleigh- Jeans tail of submm dust. ALMA 
Band 1 will allow observations in the frequency regime where PWNe dominate, and bridge 



- 14 - 



an important gap in frequency coverage, where spectral features such as power-law breaks 
occur and linear polarization observations do not suffer from significant Faraday rotation. 

4- 1.3. Radio Supernovae 

Radio supernovae occur when the blast wave of a core-collapse supernova (SN) sweeps 
through the slowly expanding wind left over from the progenitor red supergiant. Particle ac- 
celeration and magnetic field amplification lead to synchrotron radiation in a shell bounded 
by the forward and reverse shocks (Chevalier 1982). In general, free-free absorption of the 
radiation in the ionized foreground medium coupled with the expansion of the SN causes 
the radio light curve first to rise at high frequencies and subsequently at progressively lower 
frequencies while the optical depth decreases. When the optical depth has reached approxi- 
mately unity, the radio light curve peaks and decreases thereafter (e.g. Weiler et al. 2002). 
These characteristics allow estimates to be made of the density profiles of the expanding 
ejecta and the circumstellar medium and also of the mass loss of the progenitor. Resolved 
images of SNe provide information, e.g., on the structure of the shell, size, expansion ve- 
locity, age, deceleration, and magnetic field, in addition to refined estimates of the density 
profiles and the mass loss (Bartel et al. 2002). Radio observations of SNe can be regarded 
as a time-machine where the history of the mass loss of the progenitor is recorded tens of 
thousands of years before the star died. Finally, the SN images can be used to make a movie 
of the expanding shell of radio emission and to obtain a geometric estimate of the distance 
to the host galaxy (Bartel et al. 2007). 

ALMA Band 1 will allow exciting science to be done in the areas of radio light curve 
measurements, imaging of nearby SNe and, in conjunction with VLBI, imaging of more 
distant SNe. Depending on the medium, the delay between the peak of the radio light curve 
at, say, 20 cm and 1 cm, can be as long as 10 years as for instance in the case of SN 1996cr 
(Bauer et al. 2008). Absorption can also occur in the source itself. In case of SN 1986 J, 
a new component appeared in the radio spectrum and in the VLBI images about 20 years 
after the explosion and then only at or around 1 cm wavelength. The component is located 
in the projected center of the shell-like structure of the SN and may be emission from a 
very dense clump fortuitously close to that center, or possibly from a pulsar wind nebula 
in the physical center of the shell (Fig. [7| Bietenholz et al. 2004). Observations in Band 1 
minimize the absorption effect relative to observations at longer wavelengths and thus allow 
investigations of SNe at the earliest times without compromising too much on the signal to 
noise ratio of a source with a steep spectrum. ALMA with Band 1 has the sensitivity to 
measure the radio light curves of 10s to 100 SNe. In addition, ALMA with Band 1 may be 
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particularly sensitive in finding SN factories in starburst galaxies (e.g. Lonsdale et al. 2006) 
where relatively large opacities would otherwise hinder or prevent discovery. 

ALMA with Band 1 will allow high-dynamic range images of SN 1987A in the Large 
Magellanic Cloud with a resolution of about 300 FWHM beams across the area of the shell 
in 2014. This is a significant improvement over presently obtainable images (Gaensler et al. 
2007). Also, since the size of the SN increases by one Band 1-FWHM beam width per 3 
years, the expansion of the shell can be monitored accurately and in detail and will make 
this SN an important target for ALMA. 

ALMA with Band 1 could also be an important element of VLBI arrays in the northern 
as well as in the southern hemisphere. Supernova VLBI observations at 1 cm wavelength 
have provided clues about physical conditions at the earliest times after the transition from 
opaqueness to transparency, and SN VLBI with Band 1 will focus on this area of research. 

4.1.4- X-ray Binaries 

X-ray binaries (i.e., binary star systems with either a neutron star or a black hole 
accreting from a close companion) frequently show jet emission. Most of these systems are 
transients. Typically, 1-2 black hole X-ray binaries undergo a transient outburst per year, 
and the majority of such sources are at southern declinations where observations with the 
VLA are difficult, if not impossible. The rate for neutron stars is slightly larger. These 
outbursts typically last several months (although there are some which are both considerably 
longer and shorter), and over the course of the outbursts, the X-ray luminosity can change by 
as much as 7 orders of magnitude. The radio luminosities of systems seen to date correlate 
well with the hard X-ray luminosities (i.e., the luminosity above ~20 keV), albeit with 
considerable, yet poorly understood scatter. 

When the X-ray spectra become dominated by thermal X-ray emission, the radio emis- 
sion often turns off (e.g., Tananbaum et al. 1972; Fender et al. 1999), but the extent to 
which the flux turns down is still poorly constrained. This turndown is not seen in neutron 
star X-ray binaries (Migliari et al. 2004). The reduced radio emission in black hole X-ray 
binaries when they have soft X-ray spectra can be explained by models of jet production 
in which the jet power scales with the polodial component of the magnetic field of the ac- 
cretion fiow (e.g.,, Livio, Ogilvie & Pringle 1999), and may have implications for the radio 
loud/radio quiet quasar dichotomy (e.g., Meier 1999; Maccarone, Gallo & Fender 2003), and 
the still-present radio emission from neutron stars in their soft state may be indicating that 
the neutron star boundary layers play an important role in powering jets (Maccarone 2008). 
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The soft states of X-ray transients are short-lived, and during them, there may be decaying 
emission from transient radio flares launched during the state transitions. Therefore, to 
place better upper limits on the radio jets produced during the soft state, one needs a high 
sensitivity, high frequency system which is operational with a very high duty cycle. 

The radio properties of X-ray binaries with neutron star primaries are much more poorly 
understood than those of black hole X-ray binaries. This situation is partially because the 
neutron star X-ray binaries are fainter in X-rays than are the black hole X-ray binaries. 
There is, however, additionally some evidence that neutron star X-ray binaries show a steeper 
relation between X-ray luminosity and radio luminosity than do the black hole X-ray binaries, 
with oc L^x black holes, and oc for the neutron stars, a finding which 

can be explained if the neutron stars are radiatively efficient (i.e. with the X-ray luminosity 
scaling with the accretion rate), while the black holes are radiatively inefficient (i.e. with 
the X-ray luminosity scaling with the square of the accretion rate, as has been proposed by 
Narayan & Yi 1994) ~ see Koerding et al. (2006). Radio/X-ray correlations for neutron star 
X-ray binaries are, to date, based on small numbers of data points from small numbers of 
sources, and the most recent work (Tudose et al. 2009) indicates that the situation may be 
far more complex than the picture presented above. 

The faintest X-ray binaries typically show flat spectrum synchrotron radiation spectra 
from the radio through the infrared. At first blush, it may therefore appear that such sources 
arc likely to be about as easy to detect in Band 3 as in Band 1. Band 1 is essential for accom- 
plishing the best X-ray binary science possible, however, because of the strong variability 
on all timescales of these objects. Given the rarity of black hole X-ray binary transients, 
and the speed with which they evolve, one needs a system that is continuously operational. 
Additionally, including a frequency common to ALMA and the Austraha Telescope Compact 
Array would allow the possibility of continuous monitoring of X-ray binary outbursts. 

4.2. Fine Structure of Chemical Differentiation in Dcirk Cloud Cores 

Previous single-dish millimeter molecular-line observations have found that molecular 
distributions differ significantly between individual dark cloud cores. A widely accepted 
interpretation of this chemical differentiation is that there exists non-equilibrium gas-phase 
chemical evolution through ion-molecule reactions within dark cloud cores. "Young" cores 
are rich in carbon-chain molecules such as CCS and HC3N, while more evolved cores, closer 
to the protostellar formation via gravitational collapse, are rich in "latc-type" molecules such 
as NH3 and SO (Suzuki et al. 1992). Recent high-resolution mm observations, however, have 
revealed that there are even finer variations of molecular distributions down to the scale 
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of ~3000 AU, and that these fine-scale chemical fluctuations cannot be explained by the 
simple scenario for chemical evolution of cores (Takakuwa et al. 2003, Buckle et al. 2006). 
The suggestion is that there is first molecular depletion onto grain surfaces in these regions 
and then subsequent reaction and desorption of molecules back to the gas phase through 
clump-to-clump collisions or energy injection from newly formed protostars (e.g., Buckle et 
al. 2006). 

To probe the smallest length scales of chemical variations in dark cloud cores and to 
clarify the relationship among different molecular distributions (in conjunction with the 
chemical models) requires both ALMA's excellent spatial resolution and its ability to re- 
cover the large-scale structure of the cloud through observations with the ACA. As well, 
many of the most interesting molecules needed for the chemical studies within these clouds 
are readily observable across the Band 1 Receiver frequency range. These include: 



Table 1: Molecular Transitions between 34 and 49 GHz (ALMA Band 1: 31-45 GHz). 



CH3CCH 


2- 


-1 


34, 


.183414 


GHz 


HC5N 


13- 


-12 


34, 


.614386 


GHz 


so 


23- 


-22 


36, 


.202040 


GHz 


HC3N 


4- 


-3 


36, 


.392332 


GHz 


HCS+ 


1- 


-0 


42, 


.674205 


GHz 


HC5N 


17- 


-16 


45, 


.264721 


GHz 


CCS 


43- 


-32 


45, 


.379033 


GHz 


HC3N 


5- 


-4 


45, 


.490316 


GHz 


cccs 


8- 


-7 


46, 


.245621 


GHz 


C3H2 


2ii- 


-2o2 


46, 


.755621 


GHz 


c^^s 


1- 


-0 


48, 


.206956 


GHz 


CH30H 


lo- 


-Oo 


48, 


.372467 


GHz 



4.3. Mciser Science 

Maser (Microwave Amplification by Stimulated Emission of Radiation) radiation fre- 
quently occurs in regions of active star formation, where this type of emission arises from 
molecular transitions whose populations are either radiativcly or collisionally inverted. A 
photon emitted from this material will interact with other excited molecules along its path 
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causing them to emit stimulated photons with identical characteristics. This process leads 
to the creation of a highly directional beam of mono-energetic photons, which has sufficient 
intensity to be detected at very large distances. 

Masers are observed from a variety of molecular and atomic species and each of them 
serves as a signpost for a specific phenomenon, a property which renders masers powerful 



astrophysical tools (Menten 2007). More precisely, masers are formed under specific condi- 
tions, and the detection of maser emission therefore suggests that physical conditions (e.g., 
temperature, density, and molecular abundance) in the region where the maser forms lie 



within a defined range (c.f., Cohen 1995[ EUingsen 2004 and references therein). This is 



why, for example, interferometric blind and targeted surveys of maser species can lead to 



the diagnosis of interesting objects and their evolutionary phases (EUingsen et al. 2007). 



Theoretical models of masers strongly depend on these parameters as well as the geom- 
etry of the maser source. A successful model should be able to not only reproduce observa- 
tional characteristics of previously observed maser lines but also to predict new maser tran- 



sitions (e.g., the models of Sobolev et al. 1997 for Class II methanol masers and Neufeld & 
Melnick [1991 for water masers). In that respect, interferometry is essential for the successful 
search of candidate lines and being able to confirm the maser nature of the radiation. With 
its ability to provide precise images at high angular resolution, ALMA will resolve closely 
spaced maser spots and help further establish precise models of masing sources by determin- 



ing whether the detected maser signals are associated with thermal emission (Sobolev et al. 



1999). These types of information are essential for improving theoretical models. With the 



inclusion of Band 1, ALMA will cover the frequency range of 31 GHz to 950 GHz, making 
it ideal for multi-transition observations of various maser species across the mm and submm 
windows. Examples of species with observed maser radiation in the different ALMA bands 
are given in Table [2| while in Tables |3] & |4] we list SiO and methanol maser transitions that 
are observed and/or predicted at frequencies covered in the ALMA Band 1. 

Finally, maser radiation can be linearly and/or circularly polarized depending on the 
magnetic properties of the molecule. Polarimetry studies of maser radiation with interfer- 
ometers can therefore yield information on the morphology of the magnetic field threading 
the region on small scales, with the plane-of-the-sky and line-of-sight components of the 
field being probed using linear and circular polarization measurements, respectively (see for 



example Harvey-Smith et al. 2008 Vlemmings et al. 2006). Polarization data are essential 



for improving on the theory of maser polarization first introduced by Goldreich, Keeley, & 



Kwan (1973), which applies to a linear maser region, a constant magnetic field, the simplest 



energy states for a masing transition, and asymptotics limits. Observations at higher spatial 



resolution are needed to verify and improve on more realistic and extensive models (Watson 



2008). 



- 19 - 



Table 2: ALMA bands with known maser lines (Wooten 2007). 



Species 


ALMA Bands 


H2O 


B3, B5, B6, B7, B8, B9 


CH3OH 


Bl, B3, B4, B6 


SiO 


Bl, B2, B3, B4, B5, B6, B7 


HCN 


B3, B4, B6, B7, B9 



Table 3: Observed SiO maser lines in the Band 1 of ALMA (Wooten 2007) 



Transitions Frequency (MHz) 

v=0 (J= 1 ^ 0) 43423.585 

v=l (J= 1 0) 43122.079 

v=2 (J= 1 0) 42820.582 

v=3 (J= 1 ^ 0) 42519.373 

v=0 (J= 1 ^ 0) 42879.916 

v=0 (J= 1 -> 0) 42373.359 
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Table 4: Observed ( Menten ]|2007 ) and predicted (designated with a star, Cragg et al. ||2005" ) 
methanol maser lines in Band 1 



Transitions 



14(-3) -> 15(-2) E 
4(-l) ^ 3(0)E 
7(-2) ^ 8(-l)E 
6(2) ^ 5(3)A+ 
6(2) ^ 5(3)A- 
7(0) ^ 6(1)A+ 
2(0) ^ 3(1)E * 



Frequency (MHz) 



34236.8 
36169.3 
37703.7 
38293.3 
38452.7 
44069.4 
44955.8 



4.4. Magnetic Field Strengths from CCS Zeeman Measurements 

Magnetic fields are believed to play a crucial role in the star formation process. Various 
theoretical and numerical studies explain how magnetic fields can account for the support 
of clouds against self-gravity, the formation of cloud cores, the persistence of supersonic line 
widths, and the low specific angular momentum of cloud cores and stars (McKee & Ostriker, 
2007). The standard model suggests that the initial mass-to- (magnetic) fiux ratio, M/$j„jt, 
is the key parameter governing the fate of molecular cores; if M/ ^init of a core is greater than 
the critical value, the core will collapse and form stars in short time scales, but for cores with 
M/^init smaller than the critical value the process of ambipolar diffusion will take a long 
time to reduce the magnetic pressure (Mouschovias & Spitzer 1976; Shu et al. 1987). On 
the other hand, recent MHD simulations suggest that turbulence can control the formation 
of clouds and cores and in such cases the mass-to-flux ratio in the center of a collapsing core 
will be larger than that in its envelope, the opposite of the ambipolar diffusion results (Dib 
et al. 2007). Therefore, measuring the magnetic fleld strengths and the mass-to-flux ratios 
in the core and envelope separately provides a critical test for the star formation theories. 

In spite of its importance, the magnetic fleld is the most poorly measured parameter in 
the star formation process. The main problem is that magnetic flelds can be measured only 
via polarized radiation, which requires extremely high sensitivity for detections. As a result, 
the observed data on magnetic flelds is sparse compared with the observed data related to 
densities, temperatures, and kinematics in star-forming cores. The large collecting area of 
ALMA provides the best opportunity to resolve the sensitivity problem for magnetic fleld 
measurements. 
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The key to measuring the mass-to-flux ratios is the measurement of the strength of 
magnetic fields. The Zeeman effect in spectral lines is the only known method to measure 
the field strength directly. Observations of Zeeman splitting involve detecting the difference 
between left and right circular polarizations, which is generally very small in interstellar 
conditions (with the exception of masers). Successful non-maser detections of the Zeeman 
effect in molecular clouds have only been carried out with HI, OH, and CN lines because 
these species have the largest Zeeman splitting factors (~2 - 3.3 Hz///G) among all molecular 
lines (Crutcher et al. 1996, 1999; Falgarone et al. 2008). However, thermal HI and OH lines 
probe relatively low-density gas (^(H) < 10^ cm~'^) and CN is not particularly abundant 
in molecular cores. Currently, the IRAM-30 m telescope is the only instrument that has 
successfully detected the CN Zeeman effect in 7 objects. 

Our proposed project to build the ALMA Band 1 receivers provides the opportunity to 
detect Zeeman effect from CCS lines and hence greatly advance our understanding in star 
formation. The main reason is that CCS has been widely recognized as being present only 
very early in the star-forming process through chemical models (Aikawa et al. 2001, 2005) 
and observations (Suzuki et al. 1992; Lai & Crutcher 2000). Therefore the mass-to-fiux ratio 
derived from the CCS Zeeman measurements will be very close to the initial values before 
the onset of gravitational collapse. CCS also has relatively large Zeeman splitting factors (~ 
0.7 Hz//iG, Shinnaga & Yamamoto 2000) compared to most molecules. ALMA's antennas 
and site will be excellent at these "long" wavelengths, providing the stability and accuracy 
needed for such sensitive polarization work. The linearly polarized detectors on ALMA's 
antennas will also be ideally suited to measurement of Stokes V signatures from CCS. 

Using the BIMA array survey results from Lai & Crutcher (2000), Fig. [s] demonstrates 
that detections of CCS Zeeman effects can be achieved if the ALMA specifications for Band 
1 receivers are met. The detection of the Zeeman effect depends on two factors: the magnetic 
field strength and the line intensity. The two lines in Fig. [8] show the 3a detection limits for 
Stokes V spectrum with channel width of 0.024 km/s and Ihr or 10 hr integration time for 
a range of magnetic field strengths and line intensities. The channel width is chosen to have 
at least 6 channels across the FWHM of the total intensity spectrum (Stokes I). If we scale 
the line intensity from Lai & Crutcher (2000) assuming the intensity distribution is uniform 
within the 30" BIMA beam, the expected line intensity would be around 0.1-0.4 Jy for ALMA 
observations with 10" beam. Therefore, Fig. |8] shows that for the magnetic fields of 0.2-1 mC 
(typical values estimated from the application of the Chandrasehkar- Fermi method to dust 
polarimetry in dense cores), we can detect the CCS Zeeman effect with reasonable on-source 
integration time (less than 10 hr). 

In summary, building ALMA Band 1 receivers provides the opportunity to measure the 
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initial mass-to-flux ratio of molecular cores through the detection of CCS Zeeman effect, 
which cannot be done with any other instruments in the foreseeable future. The results will 
allow us to realistically test the expectations from theoretical and numerical models for the 
first time. 



4.5. Science of Jets and Outflows 

4-5.1. Jets from Young Stars 

Radio continuum emission is observed from jets and winds of young stellar objects and 
is due to the interaction of free electrons, i.e., "free-free emission". The radio images appear 
elongated and jet-like and are usually located near the base of large optical Herbig-Haro flows 
(Reipurth & Bally 2000). These regions usually have only sub-arcsecond sizes, indicating the 
youth of the emitting material and the short dynamical times involved. The emitted flux is 
usually weak, with a flat to positive spectral index with increasing frequencies (e.g., Anglada 
1995). Multi- wavelength studies of the brightest radio jets at centimeter wavelengths trace 
either earlier and stronger sources or more massive systems. The tertiary system L1551-IRS 
5, one of the most studied low-mass systems (Rodriguez et al. 1998, 2003; Lim & Takakuwa 
2006), is illustrative of the sub-arcsecond scales required (Figure [9]). 

Ground-based, interferometric studies of radio jets provide the best opportunity to re- 
solve the finest scales of the underlying source, comparable or better than optical studies 
of jets by the Hubble Space Telescope. Such fine-detail images are providing the ability to 
differentiate between theoretical ideas about the nature of these jets: i.e., the launch region, 
the collimation process, and the structure of the inner disks. Modeling efforts with the radio 
continuum emission presented in Shang et al. (2004) demonstrate one such possibility in con- 
straining theoretical parameters using present mm and centimeter interferometers (Figure [9]). 
Band 1 observations will discriminate between competing jet launch theories tied to the disk 
location of the launch point by achieving better than 0.1 arcsecond angular resolution. 

The high sensitivity of Band 1 on ALMA will also allow detection of radio emission from 
less luminous sources. ALMA will thus have the potential to discover a significant number of 
new radio jets, providing a catalog from which evolutionary changes in the physical properties 
can be deduced. As well, multi-epoch surveys will be able to follow the evolution of the freshly 
ejected material down to a few AUs from the driving sources through movies. The 31-45 GHz 
frequency range of Band 1 will show contributions to the observed emission from both the 
ionized component of the jet and the thermal emission from the dust. Together with lower 
frequency EVLA observations and detailed theoretical modelling, a complete understanding 
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of properties of the spectral energy distribution (SED) from the ionized inner regions of 
young stellar jets will be uncovered. 

4.5.2. Molecular Outflows from Young Stars 

The Sub-millimeter Array (SMA) has proven to be a successful instrument for the study 
of the youngest molecular outflows and jets from the most deeply embedded sources (e.g. 
Hirano et al. 2006; Palau et al. 2006; Lee et al. 2007a,b, 2008, 2009). The detection of 
excitation from rotational transitions of the molecule SiO up to levels J = 8-7, and CO 
up to J = 3-2 have uniquely identified a molecular high-velocity jet-hke component located 
within the outflow shell, which displays similarities to the optical forbidden line jets observed 
in T-Tauri stars (Hirano et al. 2006; Palau et al. 2006; Codella et al. 2007; Cabrit et al. 2007). 
These observations have provided a new probe of the launch of jets close to protostars and 
the mechanism of coUimation in the earliest protostellar phase. 

One unique opportunity offered by Band 1 is observation of the J = 1-0 transition of 
the SiO molecule. The SiO (1-0) transition has not yet been detected nor surveyed around 
even the brightest molecular outflows, except using single-dish telescopes (Haschick & Ho 
1990). One feature of this line that may be potentially distinct from the high- J transitions 
is that the SiO J= 1-0, like CO J= 1-0, may be tracing the outer and more diffuse gas 
that can be easily excited by shocks, located on the outskirts of outflow shells. Potential 
morphological and kinematic studies of these interesting regions where the outflows interact 
with their own prenatal clouds could be contrasted with other transitions using knowledge 
of their excitation conditions. A cooperative search along with higher transitions and with 
other envelope tracers should prove fruitful with the broad frequency coverage of ALMA. 

4.6. CO Studies of Distant Galeixies 

4.6.1. Co-Evolution of Star Formation and Active Galactic Nuclei 

Roughly half of the high-redshift objects detected in CO line emission are believed to 
host an active galactic nucleus (AGN). Although they are selected based on their AGN 
properties, optically luminous high-redshift quasars exhibit many characteristics indicative 
of ongoing star- formation, e.g., thermal emission from warm dust (Wang et al. 2008) or 
extended UV continuum emission. Indeed, galaxies with AGNs in the local Universe reveal 
a strong correlation between the mass (m) in their super massive black hole (SMBH) and 
that in their stellar bulge (measured from the stellar velocity dispersion - a; e.g. Kormendy 
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& Richstone 1995; Magorrian et al. 1998; Gebhardt et al. 2000). Such a correlation can be 
explained if the SMBH formed coevally with the stellar bulge, implying that the luminous 
quasar activity signaling the formation of a sub-arcsecond SMBH at high-redshift should be 
accompanied by starburst activity. High spatial resolution observations of CO line emission 
in high-redshift quasars can be used to infer the dynamical masses, which are found to be 
comparable to the derived molecular gas + black hole masses, meaning that their stellar 
component cannot contribute a large fraction of the total mass. There is mounting evidence 
that quasar host galaxies at redshifts, 4 < z < Q, have SMBH masses up to an order of 
magnitude larger than that expected from their bulge masses and the local relation (Walter 
et al. 2004; Riechers et al. in prep.), suggesting that the SMBH may have formed first. 
The possible time evolution of the m — a relation is of fundamental importance in studies 
of galaxy evolution, and therefore this finding needs to be put on a statistically significant 
basis. Future observations of high-redshift AGN with the Band 1 receivers on ALMA would 
allow us to address this question through the study of low-J CO line emission in galaxies 
beyond redshifts z ~ 1.5. 



4-6.2. The Molecular Gas Content of Star-Forming Galaxies at z 2 

While low-J CO line emission has only been detected in a few high-redshift objects, 
high-J CO line emission has been detected in more than sixty sources, most of which are 
classified as either submm/mm galaxies (SMGs) or far-infrared (FIR) luminous QSOs (see 
Solomon & Vanden Bout 2005 for a review). Most of these studies have been conducted 
with sensitive interferometers and single-dish facilities operating in the 3 mm band (ALMA 
Band 3), which is sensitive to higher- J CO line transitions at high redshift, as is illustrated in 
Figure |5] These lines generally trace warmer and denser gas, and so previous data may have 
led to a bias in our understanding of the molecular gas properties of high-redshift galaxies 
(e.g., Papadopoulos & Ivison 2002). The addition of 31-45 GHz Band 1 receivers on ALMA 
would open up the possibihty of studying the cold gas traced by the low-J transitions (J=2- 
1/1-0) in galaxies from moderate redshifts {z ~ 1.5) to those which existed when the Universe 
was reionized sometime before z > 6. 

Although many previous studies of CO line emission in high-redshift galaxies have fo- 
cused on those starburst galaxies and AGN undergoing episodes of extreme star-formation 
(^100 Mq yr~^), significant masses of molecular gas (> 10^° M©) have been discovered 
in more modest star-forming galaxies aX z = 1.5 — 2.0 (Daddi et al. 2008). These "BzK" 
galaxies are selected for their location in a B-z-K colour diagram (Daddi et al. 2004) and 
have star-formation rates of ~100 Mq yr^^ (Daddi et al. 2007), while their number density 
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is roughly a factor of 30 larger than that of the more extreme SMGs at similar redshifts. 
Observations of CO J— 2-1 line emission in these BzK galaxies reveal comparable masses of 
molecular gas to that of the SMGs, so their star-formation efficiencies appear lower. The 
excitation conditions of their molecular gas (temperature and density) are similar to those 
of the Milky Way (Dannerbauer et al. 2008), as indicated by the "turnover" in the CO line 
spectral energy distribution occuring at the J=3-2 transition, i.e., lower than that of the 
SMGs which typically occurs at the J=6-5 or J=5-4 transition (Weiss et al. 2005). To 
develop a full spectral energy distribution for the CO line excitation, observations of these 
galaxies in the J— 1-0 transition with the Band 1 receivers on ALMA are needed. Such ob- 
servations would also provide a more robust estimate of the total molecular gas mass, along 
with the spatial resolution needed to constrain the kinematics of the gas, as has been done 
for the SMGs (Tacconi et al. 2006). 



4.7. Substructures in Clusters of Galaxies 

Clusters of galaxies are the largest viriahzed systems in the Universe, and as such offer 
an ideal laboratory for testing how gravitational instability at large scales leads to a much 
richer range of physical phenomena at small scales. Although clusters can be studied at a 
wide range of wavelengths, a particularly powerful probe comes from using inverse Compton 
scattering of CMB photons off hot electrons in the cluster, known as the Sunyaev-Zel'dovich 
(SZ) effect (Sunyaev & Zel'dovich 1972; Rephaeli, Y, 1995; Birkinshaw,1999). This leads to 
a characteristic decrement at radio wavelengths, with an amplitude which is the integral of 
electron pressure through the cluster. 

An area of intensive current activity is to use counts of clusters found in blank sky SZ 
surveys to distinguish among models for Dark Energy. The first such detections have recently 
been announced by the SPT and ACT collaborations (Staniszewski et al. 2009; Hincks et al. 
2009), with much larger samples of clusters expected in the next few years. To turn these 
results into constraints on the equation of state of the Dark Energy it is necessary to carefully 
cahbrate mass estimators for the clusters, and also to fully characterise the selection biases in 
the samples. ALMA Band 1 could play an important role here in mapping the sub-structure 
within individual clusters to complement the single-dish data, which typically have beams 
of ~ 1', and so measure only the bulk properties of each cluster. 

However, the real power of ALMA will lie in the study of physical processes within 
the cluster for their own sake. There has been a renaissance in cluster astrophysics in the 
last several years, with an understanding of the complexity of the processes involved (see 
e.g., Markevitch et al., Giacconi et al.). The literature is now replete with phrases such as: 
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subcluster infall; cold fronts; entropy injection; cosmic ray acceleration; cold core sloshing; 
ram-pressure stripping; super-winds; jet-stirring; radio bubbles; intra-clustcr turbulence; 
magnetic draping; etc. The simple picture of virialization and hydrostatic equilibrium, which 
may apply on the average to an entire cluster (but this needs to be tested) certainly breaks 
down at smaller scales. 

ALMA Band 1 has the sensitivity to detect the SZ effect for halos from individual 
galaxies (at least for massive ellipticals), as well as for groups (postulated signal strength 
of 0.02 mJy near 30 GHz), and has angular resolution to resolve individual galaxies even 
out to 2; ~ 1 [1"= 8 kpc physical]. This means wc are entering new territory, in which the 
substructure within clusters can be probed in detail. 

Infalling sub-clusters are typically in the transonic regime, thus producing a bow shock. 
High-resolution mapping of a merger with the bow shock and its pre- and post-shocked area 
can reveal detailed microphysics of the intracluster medium (ICM, Sarazin 1988). The stand- 
off distance between a bow shock and the merger constrains the Mach number. The thickness 
of the bow shock may reflect the equilibrium time scales between ions and electrons, and 
the viscosity of the ICM, which is difficult to derive theoretically. The magnetic fields and 
surface instabilities will shape the merger during the stripping-off process. Different power 
law indices in the electron distribution, and hence in the postshock pressure, may indicate 
different particle acceleration processes in the shock front. Deriving a power spectrum from 
a pressure map tests very sensitively the underlying type of turbulence (e.g., Kolmogorov 
versus Kraichnan). A substantial level of turbulence is expected from the passage of a 
merger through the ICM (e.g.. Coma cluster core, see Schuecker et al. 2004). Finally, any 
displacement between the surface brightness peak (described by the X-ray emitting gas) and 
the pressure compared with the central cD (tracing the potential well) can indicate how 
tightly dark matter and gas are bound. Several galaxy cluster mergers have been observed 
recently with Chandra and XMM in X-rays with resolutions at the arcsecond level where 
substructures become visible (Markevitch, et al. 2000, 2002). 

Most of the above outhned physical quantities are directly related to pressure. A merger 
event will heat and compress the ICM, thus significantly increase the pressure. In the 

mm/submm bands, the Sunyaev-Zel'dovich effect (SZE, Sunyaev & Zeldovich 1972, Rephaeli 
1995, Birkinshaw 1999, Carlstrom et al. 2002) offers unique possibilities to directly probe the 
electron pressure, and it is therefore sensitive to dynamical processes. The SZE is potentially 
a more sensitive tool than X-rays, as it is directly probing the pressure ~ n T compared to 
n^y/r in the case of X-rays. It is more sensitive to high temperature regions, substructures, 
mergers and shock heated regions with pressure jumps at bow shocks. It also measures all 
the electrons and is not limited to an energy band like an X-ray observation. Furthermore, 
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it is not affected by redsfiift dimming like X-rays. In summary, tfie SZE is a promising tool 
to investigate cluster microphysics. The current highest resolution observations in SZE are 
around 13 arcseconds (Kitayama et al. 2004). 

Although the above outlined science does not require very high resolution (typically at 
the arcsecond level), relatively high sensitivity is needed to separate substructures against 
the background ICM. Clusters are typically several arcminutes in diameter. A large field 
of view is therefore desirable to avoid technical complications from mosaicing. The ALMA 
lower frequency bands arc therefore a natural choice, with Band 1 offering the largest field 
of view and the highest sensitivity and therefore a most efficient observation. At the lowest 
radio frequencies, radio point sources (which will have to be removed from high-resolution 
SZE images) are generally well characterized. 

As an illustration, two example estimates for the above science cases are given here. 
Detecting shock features - with a shock compression factor of two to three - requires resolving 
~ 100 mJy differences over ~1 arcsec. As a second example, details in the cluster gas profile 
can be detected with a ~5-10 mJy sensitivity over ~ 0.1 arcmin. 



The ALMA site is exceptionally well-suited for Band 1 observing. Even during the worst 
octile of weather, the typical optical depth through the Band 1 Receiver range is less than 
0.1. Such clear skies also benefit the other low frequency ranges, including Band 3, although 
the addition of cloud cover and water droplets in the air will make the lowest frequency 
observations more attractive. 

As shown in the previous section, the top science cases for Band 1 can stand shoulder- 
to-shoulder with the primary Level goals of ALMA. Thus, while Band 1 observations will 
clearly benefit from a larger fraction of available observing conditions, the primary motivation 
for the enhancement is not as a "poor weather" back-up receiver but rather the excellent 
science that can be attained. 



5. 



Other Considerations 



5.1. 



Weather Considerations at the ALMA Site 



5.2. 



ALMA and the EVLA 



The Expanded Very Large Array (EVLA) will also have observing capability over the 
nominal Band 1 frequency range of 31-45 GHz, through its receivers in the Ka,-band (26.5-40 
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GHz) and Q-band (40-50 GHz). Moreover, it is important to note that EVLA and ALMA 
are located at complementary latitudes, the former at +34° and the latter at —23°. With 
both telecsopes, the entire sky will be observable at high resolution at Band 1 frequencies. 
Here, we compare the relative capabilities of ALMA and EVLA over Band 1. 

Relative to ALMA, EVLA has fewer antennas (27 versus 50) but these have larger surface 
areas (25-m diameter versus 12-m), lower pointing accuracies (2-3"versus 0.6") and lower 
aperture efficiencies at Band 1 (0.34-0.39 versus 0.78). Combining these numbers (except 
pointing accuracy), the effective surface area of EVLA is a factor of 1.02-1.17 times that of 
ALMA. Adding Band 1 to the ACA antennas would likely minimize even this small difference. 
Due to its WIDAR correlator, EVLA will have the same 8 GHz maximum bandwidth as 
ALMA. The WIDAR correlator can provide users with much higher spectral resolution than 
ALMA, i.e., a maximum of below 1 Hz vs. 3.82 kHz for ALMA. Note also that EVLA will 
have continuous spectral coverage over a wider frequency range than Band 1, though there 
is some interest in pushing the upper frequency of ALMA Band 1 to 50 GHz. 

The EVLA is located at a lower altitude (2124 m) than ALMA (5000 m) and its local 
atmosphere has correspondingly larger typical precipitable water vapor. Weather statistics 
are needed for both the Plains of San Augustin and the Llano de Chajnantor to compare 
properly the availability of Band 1 observing conditions at both sites. An important factor 
that also must be quantified is the stability of the atmosphere above both sites. The ALMA 
Band 1 receiver specifications require receiver temperatures of 40-80 K over the Band, the 
same as the specification for the EVLA Ka/Q-band receivers. Assuming similar system 
temperatures for each set of receivers whenever local weather conditions are appropriate for 
Band 1 observing, the continuum sensitivities for Band 1 are similar for both observatories, 
i.e., ~1 //Jy/beam at 35 GHz after 12 hours of integration. (Variations of relative sensitivity 
across the Band 1 range have not yet been performed.) 

Given differing antenna numbers, sizes and baselines, the two observatories differ in 
various imaging metrics. (Note that these metrics ignore relative system temperatures, 
aperture efficiencies, pointing accuracies, and phase stabilities.) Comparing "single-field 
sensitivity" {ND; where D is the antenna diameter and N is the number of antennas), 
EVLA is over twice as "sensitive" as ALMA (17000 vs. 7200). Comparing "mosaic image 
sensitivity" (ND^), on the other hand, EVLA and ALMA are similar (680 vs. 600). EVLA's 
smaller number of baselines, however, yield a lower "image fidelity" (A^^) by a factor of 
> 6 (20,000 vs. 130,000). At present, EVLA has maximum baselines that are a factor of 2 
larger than ALMA's (36.4 km vs. 15-18 km), suggesting the EVLA can produce images of 
up to a factor of 2 higher resolution than ALMA can at the same frequency. EVLA's larger 
dishes, however, mean that its minimum baselines are longer than ALMA's (35-m vs. 16-m), 
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suggesting lower sensitivity to low surface brightness emission. 

In summary, for Band 1 frequencies, the EVLA and ALMA are effectively the same 
instrument at different latitudes, though each has different strengths. Weather conditions 
at both observatories suggest Band 1 observing conditions are relatively rare for EVLA but 
should be common for ALMA, but the relative fractions of available "Band 1 time" need to 
be determined. Also, the fraction of time at ALMA where no other frequency but Band 1 
can be observed needs to be determined. 



6. Summary 

The 31-45 GHz Band 1 Receiver has been considered an essential part of ALMA from 
the earhest planning days. Even through the re-baselining exercise in 2001, the importance 

of Band 1 was emphasized. With the ALMA Development Plan taking shape, we have 
undertaken an updated review of the scientific opportunity at these longer wavelengths. 
This document presents a set of compelling science cases for Band 1. The science ranges 
from nearby stars and galaxies to the re-ionization edge of the Universe. Two of the science 
cases provide additional leverage on the present ALMA Level One Science Goals and are 
seen as particularly powerful motivations for building the Band 1 Receiver: (1) detailing the 
evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires 
continuum observations out to ~ 30 GHz (~ 1cm); (2) detecting GO 3-2 spectral fine 
emission from Galaxies like the Milky Way during the era of re-ionization, 6.5 < 2; < 10 also 
requires Band 1 Receiver coverage. 
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Fig. 1 . — Spectral energy distribution plot showing the differences between three disk models 
having different maximum grain sizes. The solid curve is the model with amax = 1 cm, which 
keeps declining with roughly constant slope all the way to f cm. The two dashed curves 
are for amax = 10 yum and 1 mm. The top one, which breaks around 5 mm is the model 
with amax = 1 mm. It's interesting to note how the fluxes are very much the same for 
Omax = 1 mm or 1 cm, except precisely towards ALMA's Band 1. There is at least an order 
of magnitude difference in power at 1 cm between the maxgize = 1 mm versus the maxgize = 1 
cm disks. These models indicate that observations at the ALMA Band 1 regime are crucial 
for determining whether grain-growth to cm-sizes is indeed occurring. 
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Fig. 2. — Spectral energy distribution of TW Hya, showing the fit to the SED for an irradi- 
ated accretion disk model with a maximum particle size of 1 cm (Calvet et al. 2002). 



- 39 - 




3=^ 



i4« ^ 

[H 31 » 4« 



o 



3i.** 3SK XtO 3a » 
AlOHT ASCEMSlON {JioM] 




Fig. 3. — (Left) Image from an SPH simulation showing the surface density structure of a 
0.3 Mq disk around a 0.5 Mq star. A single dense clump has formed in the disk, at a radius 
of 75 AU and with a mass of ~ 8 Mjup. (Right) VLA 1.3 cm images toward HL Tau. The 
main image shows natural weighting with a beam of 0.11 arcseconds. The arrow indicates 
the jet axis. Upper inset: compact central peak subtracted. Lower inset: uniform weighting, 
with a beam of 0.08 arcseconds. The compact object lies to the upper right hand side. This 
condensation was also detected at 1.4 mm with the BIMA array (Welch et al. 2004). 
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Fig. 4. — VLA redshifted CO 3-2 map of the quasar J1148+5251 using the combined B- 
and C-array data sets (covering the total bandwidth, 37.5 MHz or 240 km s~^), from Walter 
et al. (2004). Contours are shown at -2, -1.4, 1.4, 2, 2.8, and 4 xcr (1 o" = 43 /iJy beam"^). 
The beam size (0.35" x 0.30") is shown in the bottom left corner; the plus sign indicates the 
SDSS position (and positional accuracy) of J1148+5251. 
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Fig. 5. — Observed frequency of ^^CO transitions as a function of redshift. The frequency 
range of the ALMA bands is also shown. 
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Fig. 6. — Three-colour image of the p Oph W photo-dissociation region (Casassus et ah 
2008). Red: MIPS 24 yum continuum Green: IRAC 8 /im continuum, dominated by the 
7.7 /im PAH Band Blue: 2MASS Ks-band image. The x— and ?/— axes show offset in RA 
and Dec from p Oph W, in degrees. The contours follow the 31 GHz emission, with levels 
at 0.067, 0.107, 0.140, 0.170, and 0.197 MJy sr-^. 
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Fig. 7. — Two-colour VLBI image of SN 1986J highlighting the emergence of a central 
component. The red colour and the contours represent the 5.0 GHz radio brightness. The 
contours are drawn at 11.3, 16. 22.690.5% of the peak brightness of 0.55 mJy/bm. The blue 
to white colours show the 15 GHz brightness of the compact, central component. The scale 
is given by the width of the picture of 9 mas. North is up and east to the left. For more 
information on the emergence of the compact source, see Bietenholz et al. (2004). 
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Fig. 8. — The expected detection limits (3 a) with integration time of 1 hr and 10 hr for a 
range of magnetic field strengths and CCS line intensity. 
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Fig. 9. — Background: The VLA+Pie Town continuum image of L1551 IRS 5 at 3.5 cm 
obtained by Rodriguez et al. (2003) in their Figure 1. The size of the beam (0.18 X 0.12"; 
P.A. = 35°) is shown in the bottom left-hand corner. Black rectangles mark the positions 
and deconvolved dimensions of the 7 mm compact protoplanetary disks. The dashed lines 
indicate the position angles of the jet cores. Inset: map of the south jet from the X-Wind 
model convolved with the beam and plotted with the same contour levels from Figure 4 of 
Shang et al. (2004). 



